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^ ' Abstract 

The problem of dual-hop transmission from a source to a destination via two parallel full-duplex relays in 
block Rayleigh fading environment is investigated. All nodes in the network are assumed to be oblivious to their 
forward channel gains; however, they have perfect information about their backward channel gains. We also assume 
a stringent decoding delay constraint of one fading block that makes the definition of ergodic (Shannon) capacity 
meaningless. The focus of this paper is on simple, efficient, and practical relaying schemes to increase the expected- 

O 

rate at the destination. For this purpose, various combinations of relaying protocols and the broadcast approach 
(multi-layer coding) are proposed. For the decode-forward (DF) relaying, the maximum finite-layer expected-rate 

> ; 

. as well as two upper-bounds on the continuous-layer expected-rate are obtained. The main feature of the proposed 
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DF scheme is that the layers being decoded at both relays are added coherently at the destination although each 
relay has no information about the number of layers being successfully decoded by the other relay. It is proved that 



the optimal coding scheme is transmitting uncorrelated signals via the relays. Next, the maximum expected-rate 
of ON/OFF based amplify-forward (AF) relaying is analytically derived. For further performance improvement, a 
hybrid decode-amplify-forward (DAF) relaying strategy, adopting the broadcast approach at the source and relays, 
is proposed and its maximum throughput and maximum finite-layer expected-rate are presented. Moreover, the 
maximum throughput and maximum expected-rate in the compress-forward (CF) relaying adopting the broadcast 
approach, using optimal quantizers and Wyner-Ziv compression at the relays, are fully derived. All theoretical 
results are illustrated by numerical simulations. As it turns out from the results, when the ratio of the relay power 
to the source power is high, the CF relaying outperforms DAF (and hence outperforms both DF and AF relaying); 
otherwise, DAF scheme is superior. 



I. Introduction 

The information theoretic aspects of wireless networks, have recently received wide attention. The 
widespread applications of wireless networks, along with many recent results in the network information 
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theory area, have motivated efficient strategies for practical applications. Fading is often used for modeling 
the wireless channels [1J. The growing demand for quality of service (QoS) and network coverage inspires 
the use of several intermediate wireless nodes to help the communication among distant nodes, which is 
referred to as relaying or multi-hopping. Many papers analyze the information theoretic and communication 
aspects of relay networks. An information theoretic view of the three-node relay channel was proposed by 
Cover and El Gamal in [|2), which was generalized in [J5] and [@] for multi-user and multi-relay networks. 
In 0, two different coding strategies were introduced. In the first strategy, originally named "cooperation" 
and later known as "decode-forward" (DF), the relay decodes the transmitted message and cooperates with 
the source to send the message in the next block. In the second strategy, "compress-forward" (CF), the 
relay compresses the received signal and sends it to the destination. Besides studying the DF and CF 
strategies, the authors in [ED-flU have studied the "amplify-forward" (AF) strategy for the Gaussian relay 
network. In AF relaying, the relay amplifies and transmits its received signal to the destination. Despite its 
simplicity, AF relaying performs well in many scenarios. El-Gamal and Zahedi jH employed AF relaying 
in the single relay channel and derived the single letter characterization of the maximum achievable rate 
using a simple linear scheme (assuming frequency division and additive white Gaussian channel). 

The problems of transmission between a disconnected source and destination via two parallel interme- 
diate nodes (the diamond channel) were analyzed in [|6l for the additive white Gaussian channels and in 
||9l for the case where the relays transmit in orthogonal frequency bands/time slots. There are also some 
asymptotic analyses on a source to destination communication via parallel relays with fading channels 
where the forward channels are known at both the transmitter and relays sides, see ifTOl and references 
therein. Diversity gains in a parallel relay network using distributed space-time codes, where channel state 
information (CSI) is only at the receivers, was presented in [fTTI . |[T2ll and references therein. 

Here, we consider the problem of maximum expected-rate in the diamond channel. A good application 
for this network is a TV broadcasting system from a satellite to cellphones through base stations where 
users with better channels might receive additional services, such as high definition TV signal lfT3ll . The 
growing adoption of broadcasting mobile TV services suggests that it has the potential to become a mass 
market application. However, the quality and success of such services are governed by guaranteeing a 
good coverage, particularly in areas that are densely populated. This paper suggests the use of relays to 
provide better coverage in such strategically important areas. The main transmitter which is a central TV 
broadcasting unit uses two parallel relays in each area with large density to improve coverage (see Fig. l]). 
According to the large number of relay pairs covering their respective areas and also the large number of 
users in each designated area, neither the main transmitter nor the relays can access the forward channel 
state information. With no delay constraint, the ergodic nature of the fading channel can be experienced 



Fig. 1. Dual-hop multicast transmission via two parallel relays. 



by sending very large transmission blocks, and the ergodic capacity is well studied [1J. According to 
the stringent delay constraint for the problem in consideration, the transmission block length is forced 
to be shorter than the dynamics of the slow fading process, though still large enough to yield a reliable 
communication. The performance of such channels are usually evaluated by outage capacity. The notion 
of capacity versus outage was introduced in [[T4l and (TJ. Shamai and Steiner [fT5l proposed a broadcast 
approach, a.k.a. multi-layer coding, for a single user block fading channel with no CSI at the transmitter, 
which maximizes the expected-rate. Since the expected-rate increases with the number of code layers 
lfT6l . they evaluated the highest expected-rate using a continuous-layer (infinite-layer) code. This idea was 
applied to a dual-hop single-user channel in [17] and a channel with two collocated cooperative users in 
lfT8l . The broadcast approach can also achieve the maximum average achievable rate in a block fading 
multiple-access channel with no CSI at the transmitters [|T9l . The optimized trade-off between the QoS 
and network coverage in a multicast network was derived in [H~3l using the broadcast approach. 

In this paper, we investigate various relaying strategies in conjunction with the broadcast approach 
(multi-layer coding) scheme for the dual-hop channel with parallel relays where neither the source (main 
transmitter) nor the relays access the forward channels. Throughout the paper, we assume that channel gains 
are fixed during two consecutive blocks of transmission. The main focus of this paper is on simple and 
efficient schemes, since the relays can not buffer multiple packets and also handle large delays. Different 
relaying strategies such as DF, AF, hybrid DF-AF (DAF), and CF are considered. In DF relaying, a 
combination of the broadcast strategy and coding is proposed, such that the common layers, decoded at 
both relays, are decoded at the destination cooperatively. Note that each relay has no information about 
the number of layers being decoded by the other relay. The destination decodes from the first layer up 
to the layer that the channel condition allows. After decoding all common layers, the layers decodable 
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at just one relay are decoded. It is proved that the optimal coding strategy is transmitting uncorrelated 
signals via the relays. Since the DF relaying in conjunction with continuous-layer coding is a seemingly 
intractable problem, the maximum finite-layer expected-rate is analyzed. Furthermore, two upper-bounds 
for the maximum continuous -layer expected-rate in DF are obtained. In the DF relaying, the relays must 
know the codebook of the source and have enough time to decode the received signal. In the networks 
without these conditions, AF relaying is considered next. Both the maximum throughput and the maximum 
expected-rate, using a space-time code permutation between the relays, are derived. In the same direction 
and for further performance improvement, at the cost of increased complexity, a hybrid DF and AF scheme 
called DAF is proposed. In DAF with broadcast strategy, each relay decode-and-forwards a portion of the 
layers and amplify-and-forwards the rest. Afterwards, a multi-layer CF relaying is presented. In the CF 
relaying, the relays do not decode their received signals; instead, compress the signals by performing the 
optimal quantization in the Wyner-Ziv sense [|20l . which means each relay quantizes its received signal 
relying on the side information from the other relay. Besides the proposed achievable expected-rates, some 
upper bounds based on the channel enhancement idea and the max-flow min-cut theorem are obtained. 
In all the proposed relaying strategies combined with the broadcast strategy, the maximum expected-rate 
increases with the number of code layers. It is numerically shown that when the ratio of the relay power 
to the source power is large, the CF relaying outperforms DAF, and hence outperforms both DF and AF; 
otherwise, DAF is the superior scheme. 

The rest of this paper is organized as follows: In Section [ill preliminaries are presented. Next, DF, AF, 
DAF, and CF relaying strategies in conjunc tion with the broadcast approach are elaborated in Sections 



to 



VII respectively. Afterwards, in Section IVIIl som e upp er bounds on the maximum expected-rate are 



obtained. Numerical results are presented in Section IVIIll Finally, Section |lX| concludes the paper. 



II. Preliminaries 

A. Notation 

Throughout the paper, we represent the expected operation by E(-), the probability of event A by Pr{^4}, 
the covariance matrix of random variables X and Y by R xy , the conditional covariance matrix of random 
variables X and Y by Hxy\w,z,-, the differential entropy function by %{•), and the mutual information 
function by Z(-; ■). The notation "In" is used for natural logarithm, and rates are expressed in nats. We 
denote / x (-) and F x (-) as the probability density function (PDF) and the cumulative density function (CDF) 
of random variable x, respectively. For every function F(x), consider F{x) = 1 — F(x) and F'(x) = dF d ^ . 
X is a vector and Q is a matrix. l nt denotes the n t x n t identity matrix. s° is the optimum solution with 
respect to the variable s. We denote the determinant, conjugation, matrix transpose, and matrix conjugate 
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Fig. 2. Network model of dual-hop transmission from a single-antenna source to a single-antenna destination via two single-antenna relays, 
the diamond channel. 

transpose operators by det, *, T , and \ respectively. U(-) and | • | represent the unit step function and the 
absolute value or modulus operator, respectively. tr(Q) denotes the trace of the matrix Q. CJ\f(Q, 1) denotes 
the complex Gaussian distribution with zero mean and unit variance. W is the Lambert VF-function, also 
called the omega function, which is the inverse function of f(W) = We w lETTl . Hi(x) is the exponential 
integral function, which is f°° ^-j-dt, x > 0. T(s,x) = t^e^dt is the upper incomplete gamma 
function, and T(s) = T(s,0). Throughout the paper, we assume that E(\Xi\ 2 ) = 1, Vi. 

B. Network Model 

Let us first restate the network model. As Fig. Q shows, the destination receives data via two parallel 
relays and there is no direct link between the source and the destination. The source transmits a signal 
X subject to the total power constraint P s , i.e., E (|X| 2 ) < P s , and the received signal at the tth relay 
is denoted by 

Y n =h re X + Z n , £ = 1,2 (1) 

The independent and identically distributed (i.i.d.) additive white Gaussian noise (AWGN) at the f th relay 
is represented by Z n ~ CA/"(0, 1), and h re ~ CA/"(0, 1) is the channel coefficient from the source to the 
tth relay. The £'th relay forwards a signal X n to the destination under the total power constraint P r , i.e., 
E (|X r J 2 ) < P r , £ = 1,2. The received signal at the destination is 

Y = hX ri + h 2 X r2 + Z, (2) 

where Z ~ CA/"(0, 1) is the i.i.d. AWGN and h e ~ CAf(0, 1) is the channel coefficient from the f th relay 
to the destination. All h re and hg are assumed to be constant during two consecutive transmission blocks. 
Obviously, channel gains ai = \hi\ 2 and a re = \h re \ 2 have exponential distribution. 

Note that the transmitter as well as both relays and the receiver are equipped with one antenna. We 
assume that the relays operate in a full-duplex mode and they are not capable of buffering data over 
multiple coding blocks or rescheduling tasks. Since there is no link between the relays, the half-duplex 
mode is a direct result of the full-duplex mode with frequency or time division. 



6 



C. Definitions 

In the following, the performance metrics which are widely used throughout the paper are defined. The 
expected-rate IZf is the average achievable rate when a multi-layer code is transmitted, i.e., the statistical 
expectation of the achievable rate. The maximum expected-rate, namely 7Z m , is the maximum of the 
expected-rate over all transmit covariance matrices at the relays, transmission rates in each layer, and all 
power distributions of the layers. Mathematically, 

K 

7Vf= max ^2Vi(Ri)Ri, (3) 
tr(Qi)<Pi i=1 

where Ri, Qj, and Vi are the transmission rate, transmit covariance matrix at the relays, and probability 
of successful decoding in the z'th layer, respectively. 

If a continuum of code layers are transmitted, the maximum continuous-layer (infinite-layer) expected- 
rate, namely 7Z m , is given by maximizing the continuous-layer expected-rate over the layers' power 
distribution. 

When a single-layer code is transmitted at the source and the relays, the average achievable rate is called 
the throughput, namely 7Z S . The maximum throughput, namely lZ m , is the maximum of the throughput 
over all transmit covariance matrices at the relays Q, and transmission rates R. Mathematically, 

71™ = max V(R)R. (4) 

R,Q 
tr(Q)<P 

III. Decode-Forward Relays 

In order to enhance the lucidity of this section, single-layer coding is studied first. The idea is then 
extended to multi-layer coding. Since the continuous-layer expected-rate of th is scheme is a seemingly 



intractable problem, a finite-layer coding scenario is analyzed in Section 



III-Bl. 



A. Maximum Throughput 

In single-layer coding, a signal X = 7X1 with power P s and rate R = ln(l + P s s) is transmitted, where 
7 2 = P s . The tth relay decodes and forwards the received signal in case a re > s. If a n < s, then a re is 
replaced by zero. The coding scheme at the relays is a distributed block space-time code in the Alamouti 
code sense Il22l . At time t, the first relay sends aX\(t) while the other relay sends /3Xx(t+ 1). To satisfy 
the relays power constraint, it is required that a 2 = (5 2 = P r . At time t+1, the first and the second relays 
send — aX{{t + 1) and f3X*(t), respectively. The relay with a Tl < s simply sends nothing. Applying the 
Alamouti decoding procedure and decomposing into two parallel channels, the throughput is given by 

P, 



71 



D,s 



Pr {a n > s} Pr {a T2 > s} Pr <^ a\ + 02 > s- 
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Pr {a ri > s} Pr {a r . 2 < s} Pr { a-y > Sy } + 



Pr {a ri < s} Pr {a r2 > s} Pr < a 2 > s 



ln(l + P a s). 



(5) 



The first term in the right hand side of Eq. <J3|) represents the case of decoding the signal at both relays 
and the destination. The second and third terms represent the probability of decoding the signal at only 
one relay and the destination. Substituting the channel gain CDFs in ©, the throughput is given by 



y 



Ud ' s = (^e- s -e^ + 2^e- s (* +1 )ln(l + F sS ). 



(6) 



Theorem [1] proves the optimality of the above scheme and presents the maximum throughput of the 
channel. 



Theorem 1 In the proposed single-layer DF, the maximum throughput is achieved by sending uncorrelated 
signals on the relays, the maximum throughput is given by 



, Ds — max — se 

U ' S 0<s<s t \P 7 



(^-se- s -e- s + 2y-<^ +1 ) ln(l + F s s), 



(7) 



where s t = min |2^, 1.212 j. 
Proof: 



A 



Consider Q = P, 



1 P 
P 1 



as the relays transmit covariance matrix. Therefore, E (X ri X* 2 ) = pP r . In 



the following, we shall show that p° = 0. Let us define F(s) as follows 

F(s) = Pr {a ri > s} Pr {a r2 > s} Pr ja > s^- 



Pr {a ri > s} Pr {a r , 2 < s} Pr ^ a x > s— f + 



Pr {a ri < s} Pr {a T . 2 > s} Pr < a 2 > s 



(8) 



where a = -^hQh) and h 
is 



hi hi 



The maximum throughput of the diamond channel in general form 



Tl% s = max F(s) ln(l + P s s) 



(9) 



The only term in F(s) which depends on p is Pr ( a > Sjr )• Since Q is non-negative definite, one can 



write it as Q = UDTJt, where D = P r 



1 + p 
1- p 



is non-negative diagonal and U = 4j 



1 1 



is unitary. Since hi and h 2 are independent complex Gaussian random variables, each with independent 



zero-mean and equal variance real and imaginary parts, the distribution of hXJ is the same as that of h 
G31 Thus, 



Pr | a > s^j = Pr |£q# > P s s} 

= Pr | (Ku) D (hU^ > P s s| 



Pr ^KDh) > P s s| . (10) 



The last expression in Eq. (| 1 Oh corresponds to the complementary CDF in MISO channels. Jorswieck 
and Boch ll24ll proved that in an uncorrected MISO channel with no CSI at the transmitter, but perfect CSI 
at the receiver, for every transmission rate, the optimal transmit strategy minimizing the outage probability 
is to use a fraction of all available transmit antennas and perform equal power allocation with uncorrected 
signals. Therefore, the solution of max p tr ( D )> 2 p r Pr jln ^1 + fo'D/iM > hi(l + P s s)| is p = or p = 1. 



Defining 



S ' = -( 2W -(2^) +1 )f SS2 ' 512 4 <U) 



where WLi (•) is the -1 branch of the Lambert W-function, one can show that if s < s c , then 

F p=0 (s) > F p=1 (s). (12) 

In the remainder of the proof, we shall show that in case p = 1, s° < s c . Then, as Vs < s c , F p=0 (s°) > 
F p= i(s°), it implies p° = 0, i.e., the optimum correlation coefficient between the relay signals maximizing 
the throughput of DF diamond channel is zero. 

Assume that s° maximizes 1Z(s) = F p= i(s) In (1 + P s s). Hence, lZ'(s°) = 0. Defining f p= i(s) = 
-F p=1 (s), we get 

K'{s) = F p=1 (s)—^— - f p=1 {s) In (1 + P s s) . (13) 

1 i s S 

Let us define g (s, P s ) = In (1 + P s s) p » and r(s) = • As such, we get 

^'(s)>0 iff r(s)> g(s,P s ), 

TZ'(s) = iff r(s) = g (s, P s ) , (14) 
K'{s) < iff r(s) < g(s,P s ) . 



Ps s\ . P. 



Noting F p= x(s) = ( e~ s + 2 (1 - e~ s ) e~^^\ e'^^) 8 , we have 



It can be shown that as far as s > s t = min |2jf , 1.212|, we have 

r(s) < s, Vs > St. 



The derivative of g (s, P s ) over P s is 

dg(s,P s ) sP s -\n(l+sP s ) 



if 1 



(16) 



(17) 



ap, p2 p- ' 

Therefore, g(s,P s ) is a monotonically increasing function of P s and its minimum is in P s = 0. As a 
result, 



l + P S s 



o (s, P s ) > lim In (1 + P s s) p ° = s. 



(18) 



Comparing Eq. dl6J), Eq. (jl8|), r(0) 



2P r 



Applying Eq 



to Eq. (|14t) gives 



> and g (0, P s ) = yields 

r(s) > c/(s,P s ) s = 0, 
r(s) < g(s,P s ) s > s t . 



TZ'(s) > s = 0, 
7e'0) < s > st. 



(19) 



(20) 



As TZ(s) is a continuous function, according to Eq. (|20r) . < s° < s t . Noting s 4 < s c , Eq. (|12h yields 
P p= n(s ) > F p=1 (s°) and as a result, p° = and a = a x + a 2 . Substituting the channel gain CDFs in 
Eq. d8|), the maximum throughput of the DF diamond channel is given by Eq. O), which is achievable 
by applying the aforementioned distributed space-time code. 



B. Maximum Finite-Layer Expected-Rate 

For the lucidity of this section, the encoding and decoding procedures are presented sparately. 
1) Encoding Procedure: The transmitter sends a f^-layer code X = X)i=i7i-^Q to the relays, where 
7 4 2 represents the power allocated to the i'th layer with rate 

The relays start decoding the received signal from the first layer up to the layer that their backward 
channel conditions allow. Then, the relays re-encode and forward the decoded layers to the destination. 
To design the transmission strategy, we first state Theorem Q. 



Theorem 2 In multi-layer DF, if the layers ' power distribution in the first relay is equal to that of the 
second relay, the relay signals must be uncorrelated in order to achieve the maximum expected-rate. 
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Proof: 

Analogous to the proof of Theorem Q, let us define 



+ F ari ( Si )F ar2 (8i)Vi,i + F ari ( Si )F ar2 ( Si )Vi,2, 



(22) 



where Vi } i s 2, 7\i, and 7\2 are the probability of decoding the z'th layer at the destination when both 
relays, only the first relay, and only the second relay decode the signal, respectively. The expected-rate 
in the z'th layer can be written as 



n i (s)=VAn 1 + 



The only term in Eq. (|22l) which depends on the transmit strategy at the relays is 7^ 12. We denote Q 
as the transmit covariance matrix of the relays in the z'th layer. So that, 



(23) 



Pr < 1 



EL 



> Ri 



(24) 



Analogous to the proof of Theorem ll[ by decomposing Qj and J2f=i+i Qj'» and noting the fact that 
multiplying h by any unitary matrix does not change the distribution of h, we get 



V, 



i,l,2 



Pr < 



1 + 



Pih 


1 + pi 
I - Pi 




l + hh 


1 + Pi 
1 - pi 





> Rr 



> . 



(25) 



It can be shown that the optimum solutions for p and p to minimize 7-^1 2 in Eq. (1251) is either p t = fa = 
or pi = f>i = 1 11251 . We shall now show that the optimum solution is p° — p\ — 0. Towards this, we 
follow the same general outline to the proof of Theorem Q. 
Let us define the following functions, 

(l + LsAh+a + PAsi 



9 (*^i> Pi 1 



In 1 



Pi$i 

1 + Lsj 



nsi 



Vi_ 



(26) 
(27) 



One can simply show that Eqs. 
replaced by Sj. 



fl) 



14J) and 06J) still hold by redefining the functions as above, and with s 



Defining P — 1 f* , from Eq. (18) and noting IjSj > 0, we have 



g(si,Pi,Ii) 



i + 5 ln l 1 + TT77J 
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> 



( 1 I PjSj \ 

lA" V 1 + W 



In 1 + Ps, 



(i+f°») 
p 



> Si- 



(28) 



Therefore, Eqs. (|19|) and (|20|) still hold with the above functions, and then, < s° < s t . Noting s t < s c 
results because as pointed out earlier P i p=0 (s°) > Vi jP =i(s°). 



With respect to Theorem |2j the following transmission scheme is proposed. Assume that the first and 
the second relays decode M and N layers out of the whole K transmitted layers, respectively, according 
to their corresponding backward channel. As the relays do not know the channel of the other relay, and 
hence, do not know the layers' power distribution in the other relay, its code construction is based on 
a similar power distribution assumption for the other relay. Theorem Q demonstrates that uncorrected 
signals must be transmitted over the relays. For this purpose, the following scheme is proposed. At time 
t, the first relay sends Ylf=i a iXi(t) while the other relay sends Yli=i P%Xi(t +1). At time t+1, the first 
and the second relays send Ylf=\ —otiX*(t + 1) and J2?=i PiX£(t), respectively. Note that J2iLi a i = Pr, 
cti = for i = M + 1, K and J2?=i = P r , A = for i = N + 1, K. 

The received signal at the destination is 

Y(t) = h J2*Li oaXiit) + h 2 Zf =1 PiXt(t + 1) + Z(t), 

r(t+i)=-^iE£i«^+i)+^E£iA^;W+^(t+i)- 



(29) 



One may express a matrix representation for Eq. (|29|) as 



Y(t) 
-Y*(t+1] 



K 



X t {t) 

Xi(t+r 



+ 



Z(t) 

-z*(t+r 



(30) 



h\OLi h 2 p- 
-h* 2 f3i h\a 

2) Decoding procedure: The destination starts decoding the code layers in order, from the first layer 
up to the highest layer that is decodable. To decode the z'th layer, after decoding the first % — 1 layers, the 

h*ai -h 2 P 



channels are separated into two parallel channels by multiplying both sides of Eq. (|30l ) by 
Therefore, 

Xi(t) 
Xi(t+1) 

h\oti —h 2 Pi hictj h 2 /3j Xj(t) Z(t) 
h* 2 (3i hxoii -h* 2 (3j h\aj X^t+1) Z(t+1] 



h* 2 Pi h x a.i 



' Y(t) ' 




Y(t+1)_ 





a l af + a 2 pf 

a-iaj + a 2 f3f 



+ 



K 

E 

i=i+i 



(31) 



Z(t) and Z(t + 1) are two independent i.i.d AWGN, each with power aiaf + a 2 j3f. 
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The interference power caused by upper layers while decoding the z'th layer is 

K 

j=i+i 

K 

j=i+i 

Thus, the probability that the i'th layer can be successfully decoded at the destination is 

*>, = Pr/ ^^±^1 > p ). (33 ) 

1 1 + E i= i+i (oiog + a 2/3|) i + Ej=if 1 tj* J 

Hence, the achievable expected-rate using this scheme can be written as 



K 



Knj = J>m ( 1 + TT^ZT 1 2 ) • (34) 
To summarize, we have shown the following. 

Theorem 3 In the diamond channel, the above result implies that the following expected-rate is achiev- 
able. 

K < -2, 

(35) 



TVB f = max YVi\n[l+ ^fi 



Si 



with V t = Pr | 1+E ^^Xl 2 /^) " i+SjS.^ J- ™ e ™jmtotiton to E£i7, 2 = Ps, 

Yli=i a i = J2i=i Pi = Pr, where a-i and are zero for the layers which are not decoded at the relays. 
Note that c^s and fas are optimized separately. 

Remark 1 One important feature of the proposed scheme is that the layers being decoded at both relays 
are added coherently at the destination although each relay has no information about the number of 
layers being successfully decoded by the other relay. 

IV. Amplify-Forward Relays 

A simple but efficient relaying solution for the diamond channel is to amplify and forward the received 
signals. In order for the destination to coherently decode the signals, it employs a distributed space-time 
code permutation along with the threshold-based ON/OFF power scheme, which has been shown that 
improves the performance of AF relaying [11]. According to the ON/OFF concept, any relay whose 
backward channel gain is less than a pre-determined threshold, namely a t h, is silent. In this scheme, 
the relays transmit the signals to the destination in two consecutive time slots. In time slot t, the first 
(resp. second) relay transmits c\Y n (t) (resp. c 2 Y T2 {t^ 1)). In time slot t+1, the first (resp. second) relay 
transmits —c 1 Y* l {t-\- 1) (resp. c 2 Y* 2 (t)) with the backward channel phase compensation [fTTTl . To satisfy 
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the relays' power constraint, it is required that q = \j u ^ e p a ^ Pr , £ = 1,2, where W(-) is the unit step 
function. At the destination, the channels are parallelized using the Alamouti decoding procedure ll22l . 
The received signal at the destination is 



(36) 



y(t)=ci/n(t)y ri (t) + c 2 h 2 {t)Y r2 {t + 1) + z{t), 

Y(t + 1) = - dh^Y^t + l) + c 2 h 2 (t)Y r * 2 (t) + Z(t + l). 
As the destination accesses the backward channels, after compensating the phases of h ri and h r2 into h* 
and h* 2 in time slot t + 1, we get 

h ri hiCi h r , 2 h 2 c 2 X(t) 
-K 2 h* 2 c 2 h* r h\ Cl X(t + 1) 

c x h x Z rx {t) + c 2 h 2 Z r2 (t + 1) + Z(t) 
Ci/i*Z ri (t + 1) - c 2 /i*Z r2 (t) - + 1) 
t 



Y(t) 




-Y*(t + 1) 





+ 



(37) 



Multiplying 



h ri h\Ci h r2 h 2 c 2 



-h* r2 h* 2 c 2 h* ri h\ci 
source-destination instantaneous mutual information is 



l(X;Y) = ln 1 + 



to both sides of Eq. (1371) . two channels are parallelized, and the 



(l/i n /ll 


2 c\ 


+ 


h r2 


h 2 


2 4)p s \ 


1 + 


\hx 


2/2 


+ 


h 2 


2 4 J 



which is equivalent to a point-to-point channel with the following channel gain, 



a AF,2 



(38) 



(39) 



1 + i^t+i ai + ^7fn a 2 

If one relay is silent and only one relay transmits, let say the £'th relay, by replacing zero instead of one 
of the channel gains into Eq. (I39l) . we get 

atP r 



O-AFA 



1 + a n P s + a t P r 

The expected value of the optimum ON/OFF threshold in which <iaf,2 > &af,i is given by 



ath 



1 + P s + P r 



(40) 



(41) 



Proposition [l] yields the maximum achievable throughput in this method. 
Proposition 1 The maximum achievable throughput in the above AF scheme is specified by 



K 



A, s -max e i+p.+*t [e i+ft+*v F &A (s) 



+2(1- e" t+^w F aAF : (a) m(l + P s s) , 



(42) 



where F &AF2 {-) and F &AF1 {-) are the CDFs of claf,2 and claf,i from Eqs. (I39h and (I40I) . respectively. 
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The maximum continuous-layer expected-rate of the above AF relaying is presented in Theorem ^ 



Theorem 4 The maximum achievable expected-rate in the above AF relaying is given by 

Pr x 7_ / (2 f(s) 



Pr 



Ti^=e~ ^+Ps+Pr\2 - e -i+Ps+Prj J F(s) ( - + : ^77ff ) ds, (43) 

so 

with 

= 1 , (44) 

2e 1 + p «+' p r — 1 

/(^) = n^)- (45) 

77ze integration limits are the solutions to F(sq) = sq(1 + P s so)/(so) and F(si) = sif(si), respectively. 
Proof: 

The maximum achievable expected-rate at the destination can be expressed by 

K2 c =2e- ath (1 - e~ ath ) K? + e~ 2ath TZ^ 

where VJ^ and are the maximum expected-rates when only one relay is active and both relays are 
active, respectively. As showed in lfT5l . []26l . 1Z™ and 1Z™ are given by 

, , , -sl'(s) 



TV? 


= ma*/ 


I(s) Jq 











1 + sl{s) 



ds, 



Substituting the above equations in Eq. (1461) . we get 



oo 

1 - -4 -F. c 

\ 2- e~ a *h 



3.AFA \ S ) 



-/■:,„») , aJ/ r ( f\ d.s. (48) 



Defining 



a 2(1 — p- a *h ) p-ath 

f M = L.-., W') + J— « (49) 



2 _ e -ath W - 1V ' 2 - e- a *" 
the maximum expected-rate of the proposed AF scheme is found by 



r°°— -s/'Cs) 

Kl=maxe- ath (2-e- ath ) F(s) )r/-ds. (50) 

' /(.) v J J 1 + s/(s 
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Substituting a th by l+ p r +P and maximizing over I(s) by solving the corresponding Euler equation 11271 , 
we come up with the maximum expected-rate as 

Sl 

^= e ™(2 - e - 1+^^(5)0 + ^^5, (51) 

so 

where s and s% are the solutions to F(s ) = s (l + P s so)f(s ) and F(si) = s\f(si), respectively. 

■ 

Remark 2 In the above results, the power constraint P r has been applied only to the time slots when the 

relays are ON. Alternatively, one can assume that the relays have the ability to save their power while 

working in the OFF state and consume it in the ON state. In this case, all the above calculations hold 

p r 

except for the integration limit s which is now the solution to F(sq) = sq(1 + e l + p *+ p r P s s )/(sq). 

V. Hybrid Decode-Amplify-Forward Relays 

In this section, we propose a DAF relaying strategy which takes advantage of amplifying the layers 
that could not be decoded at the relays in the DF scheme. Specifically, each relay tries to decode as many 
layers as possible and forward them by spending a portion of its power budget. The remaining power is 
dedicated to amplifying and forwarding the rest of the layers. 

In order to enhance the lucidity of this section, single-layer coding is studied first. The idea is then 
extended to multi-layer coding. As the continuous-layer expected-rate of this scheme is a seemingly 
intractable problem, a finite-layer coding scenario is analyzed. 



A. Maximum Throughput 

A single-layer code X = yXi with power P s , i.e., 7 2 = P s , and rate R = ln(l + P s s) is transmitted. 
If a ri > s, then the £'th relay decodes the signal and forwards it, otherwise, it amplifies and forwards 
the received signal to the destination. In time slot t, the first (resp. second) relay transmits X ri (t) (resp. 
X T2 (t + 1)). In time slot t + 1, the first (resp. second) relay transmits — X* (t + 1) (resp. X* 2 (t)) with the 
backward channel phase compensation. There are three possibilities: 

1) a ri > s and a r2 > s: both relays decode the signal. In this case DAF is simplified to DF in 



Section 



III. 



2) a ri < s and a r2 < s: none of the relays decodes the signal. This case is simplified to AF in 



Section 



IV 



3) a ri > s, a r . 2 < s or a ri < s, a r2 > s: only one relay decodes the signal. 
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In the third case, without loss of generality, assume that the first relay decodes the signal and the second 
relay does not decode it, i.e, a n > s,a r2 < s. Hence, X n (t) = aXi(t) and X r2 (t) = c 2 Y r2 (t) = 
c 2 (h^jX^t) + Z r2 (t)), where a 2 = P r and c 2 = 



a r „P s +l 



At the destination, we have 



Y{t) = haXtit) + h 2 c 2 h r2l X l {t + 1) 

+h 2 c 2 Z r2 (t+ 1) + Z(t), 
Y(t + 1) = -h ia X*(t + 1) + h 2 c 2 h* nl Xl{t) 

+h 2 c 2 z; 2 {t) + z{t + i). 

After compensating the phase of h r2 into h* in time slot t + 1, we get 



(52) 



Y(t) 






K 2 h 2 c 2 ^ 




X(t) 


-Y*(t + 1) 






c 2 7 h\a 




X(t + 1) 






+ 


c 2 h 2 Z r2 (t + 1) 


+ Z(t) 






—c 2 h\Z T2 (t) - 


- Z*(t+1) 



(53) 



Multiplying 



h x a h r2 h 2 c 2 'y 
-h* h* 2 c 2 -i h\a 



t 



to both sides of Eq. (|53|) . two channels are parallelized and the 



source-destination instantaneous mutual information is 

1{X;Y) =ln (1 + 



(\hi 


2 a 2 + 


h r2 h 2 


2 c 2 2l 2 ) 


1 + 


h 2 


2 r 2 
c 2 



(54) 

A comparison of this method and the DF scheme reveals that if a r2 > ^-a\, then DAF outperforms DF, 
otherwise, we switch to DF, that is the second relay becomes silent. Since the relays do not know the 
value of a\, they use its expected value. As a result, the amplification coefficient of DAF can be written 

as ce = 



U(a ri 



.It can be shown that the maximum throughput of this scheme is given by the 

dj- ^ is "T* J- *~ 

following proposition. 



Proposition 2 The maximum throughput of the proposed hybrid decode-amplify-forward relaying is given 
by 

-2e s -%-l)e-< 2+ %) 



K 



DA,s~ 



max 

s 



2e s + 2e~% + s^- 



+ (e' ath F^ 2 (s) + (l-e- a ^)F w (s))e-^(l-e- s y 



+ 2e-( s+ ^) {l-e- s )F, 



a DAF 



Ps 
Pr 



ln(l + P s s), 



(55) 



, ai+a r „ P s (ai+a.2) 

where a DAF = l+a ^ Ps+a2 p r , a th 



l + Ps+P, 



-, and qaf,2 and clafi are from 1391 and \40\ respectively. 
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B. Maximum Finite-Layer Expected-Rate 

Since continuous-layer coding for DAF relaying can not be directly solved by variations methods, 
we choose a finite-layer code and proceed as follows. In the finite-layer broadcast approach, the source 
transmits a K layer code X = J2?=i Ji^i t0 tne relays, where jf represents the power allocated to the 
i'th layer with rate 



Ri = In 1 + 



1 + EL+1 7?Si 



(56) 



<?=«+ i i j ' 

Each relay decodes its received signal from the first layer up to the layer that its backward channel 
conditions allow and forwards them to the destination. Afterwards, each relay amplifies and forwards the 
remaining undecoded layers. 

Suppose that the first and second relays allocate portions £P r and (P r of their power to the decoded 
layers, respectively. Also, assume that the first and second relays respectively decode M and N layers 
out of the K transmitted layers. Without loss of generality, assume M > N . Denote by af (resp. (5f) 
the power allocated to the z'th layer at the first (resp. second) relay. The amplifying coefficients are 



~J Pr — o — for the first relay and c 2 = W — „ — 5 — for the second relay. Let us define 



on 



A 



/i ri ci7i for % = M +1, K and = h r2 C2ji for i — N+l, K, The coding scheme is as follows. 
At time t, the first relay sends Ylf=x a iXi(t) while the other relay sends J2f=i PiX%{t + 1). At time t + 1, 
the first and the second relays send J2?=i —0(iX*(t + 1) and E;=i w i tn compensating the phases 

of h ri and h r2 into h* x and h* , respectively. 
The received signal at the destination is 

y{t)=hiY,ti ^x^+hzl, PiMt + 1) 

+h 1 c 1 Z ri (t) + h 2 c 2 Z r . 2 (t + 1) +Z(t), 

Y(t + i) = -hZf =1 a4X*(t + i) + h 2 zf =1 0*X*(t) 
-h lCl Z* ri {t + 1) + h 2 c 2 Z* 2 (t) +Z(t + 1). 



(57) 



One may express a matrix representation for Eq. (|57l) as 



Y(t) 
-Y*(t + 1) 



K 



E 



8=1 



Xi(t+1) 



+ 



h x OLi h 2 /3i 
-h* 2 (3i h\Ui 

h lCl Z ri (t) + h 2 c 2 Z r2 (t + 1) +Z{t) 
h\c x Z T ^{t + 1) - h* 2 c 2 Z r2 (t) -Z*(t + 1) 

The destination starts decoding the code layers in order, from the first layer up to the highest layer that 
is decodable. To decode the i'th layer, after decoding the first i — 1 layers, the channels are separated into 



(58) 
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two parallel channels by multiplying both sides of Eq. (1581) by 



hat h 2 (3i 



t 



. Therefore, 



" Y(t) ' 




Y(t + 1)_ 





ai«f + a 2 /3f 

axaf + a 2 (3f 



K 

E 



^•(t + l) 



X»(t + 1 

Z(t) 
Z(t+1) 

2 



+ 



+ 



(59) 



Z(t) and Z(i + 1) are two independent i.i.d. AWGN, each with power (a\oq + a 2 p~) (1 + (i\c\ + a 2 c 2 ) 
The interference power caused by upper layers while decoding the z'th layer is 

K 

k = 22 (i a l a i a j + a 20iPj) 2 + a l a 2 (<*ifij - a jPi) 2 ) 



j=i+l 



K 



= (aiOij + a 2 (5f) 2_, { a i a j + a 20f) ■ 
j=i+i 

Thus, the probability that the i'th layer can be correctly decoded at the destination is 



(60) 



->- 



[l+a 1 c 2 1 + a 2 c 2 2 +J2 j=t+1 (a 1 a 2 j + a 2 /3f) l+J2 j=l +i 7j s i) 
Hence, the expected-rate at the destination using this scheme can be written as 

i=0 V 

To summarize, we have shown the following. 



(61) 



(62) 



Theorem 5 The maximum achievable expected-rate in the proposed DAF relaying is given by 

K 



n 



max In ( 1 + ^ 

Z=0 \ ^ Z^j=?.+1 /j 



where 



DAJ 



V l= Pr- 



aia? + a 2 /3f 



->- 



l+aic2 + a 2 c|+E i=i+ i(ai^ 2 + a 2 /3f) l+£ 7?* 



(63) 



(64) 



and a, 



7i, ? = M + 1, K, and 



CPr 



«r 2 Ei = iV + l7» 2 + l 



7i, i = iV + 1, 77ze 



power constraints are Yl^ill = Ps> X^=i a l = ^Pr, and Yli=i (3f = (P r . Similar to DF scenario, 
(cxi, oi 2 , . . . , ctM, and (f3i, (3 2 , . . . , /3n, C) are optimized separately. 
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VI. Compress-Forward Relays 

In CF relaying, the relays quantize their received signals using an optimal Gaussian quantizer with 
minimum mean-square error (MSE) criterion [|28"1 , and then forward the quantized signals. With respect 
to the correlation between the relays signals, Wyner-Ziv compression method ||20l is applied. In this 
scheme, the relays do not decode the signal and hence, the latency and complexity is lower in comparison 
with DF and DAF. Also, the relays do not need to access the source codebook; however, the source-relay 
channel gains must be available at the destination. 

Denote by q ri and q r2 the quantized signals at the first and second relays, respectively. One can write 
the following equations on q re , I = 1,2, 

Y rt = q ri + n rv (65) 

and 

q n = 9iY n + n n , (66) 

where n re ~ CJ\f(0, D e ) and h re ~ CA/"(0, 9(.Di) are the equivalent quantization noises independent of q n , 
®i = 1 — i i Dl p > an d Dp is the quantizer distortion at the tth relay ||29l . 

If the destination decodes q n and q r2 , and the transmission rate is below X(X; q n ,q r2 ), the signal is 
successfully decodable. For simplicity, let us assume that the optimum value of the quantizer distortion 
D° t and the optimum value of the relays rate R° are selected independent of the source-relays channel 
gains. Hence, with respect to the network symmetry, D° = D% and R° = R° , and therefore, they are 
simply denoted by D and R r , respectively. 

To decoded the quantized signals at the destination, based on the multiple-access capacity region ll28l 
in the second-hop, the following inequalities must be satisfied, 

R r < l(X n ;Y\X r2 ) = In {a x P r + 1) , 
R r < l{X r2 ;Y\X ri ) = In (a 2 P r + 1) , 

2R r <T(X ri ,X r2 ;Y) = ln((ai + a 2 )P r + 1) . (67) 

For lossless compression of the quantized signals, based on the Wyner-ziv rate region [1201 . we have 
the following inequalities, 

R r >X(q ri ;Y ri \q r2 ), (68) 
R r > I(q r2 ]Y r2 \q ri ), (69) 
2R r >l(q n ,q r2 ;Y n ,Y r2 ). (70) 
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In the problem in consideration, Eq. (1701) is 



■£\Qri i Qr2 ) ^ri i ^2 ) ^ 



hi 



det Ry x y 2 

det RYiY 2 |q ri ,g r2 

(g ri + a f2 )P s + 1 
D 2 



(71) 



In order to derive a closed form expression for Eqs. (1681) and (1691) . let us first estate the following 
lemmas. 



Lemma 1 The mutual information between the source signal and the relays quantized signals is given 
by 

Z(q n ,q r2 ;X) = m(l + a C FP s ) , (72) 

1 ^ e 2 +i 0i 1 fli+i e 2 
Proof: The mutual information between the source signal and the relays quantized signals can be 

expressed by 

X( qr ,, qr . 2 ;X)=ln( V (74) 



det Rg ri q r2 |X 



where 



and 



Thus, 



det R qriqr2 = 9lela ri P s + 9\9 2 Da rx P s + 9\9\a n P s + 9\0\ 



2 2 L> + 9 x 9\a^P s D + 0^1) + ^6> 2j D 2 , (75) 



det R qriqr2 \x = e\6l + 0?0 2 £> + 9 X 9\D + fl^D*. (76) 



9i9 2 a ri +9iDa ri +9i9 2 a r2 +9 2 a r2 D 
Inll+K, + 



0i0 2 + + 9 2 D + D 2 

fli0 2 + 9 2 D 
9 X 9 2 + 9 X D + 9 2 D + D 2 



+ Qr 2 — . „ ^ . „ ^ . I Pi 



In I 1 + I e2D+D 2 + 6>iD+p2 I 

1 " l " 0102+6*1 1 " l " 0102+02 . 

Mi+I rrrni + TTMl lfil. C77) 

02 + 1 01 01 + 1 02 
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Equation (1731) together with Eq. (|77l) results 



Lemma 2 In the problem of interest, we have 



X(q n ; Y ri | q r% )= In (1 + a CF P s ) + In 



fll + fl) (fl 2 + £>) 

D(l + a r2 P s ) 



(78) 



Proof: 



X(q n ; y n | q r2 ) = X(q ri ; X, Y ri \ q r2 ) - X(q n ;X\Y ri , q r2 ) 

— Z(q ri ;X, Y ri \q r2 ) 

— ^{ ( lr 1 , Qr 2 ! X, Yri ) ~ ^{Qr 2 X, Y ri ) 

= X(q ri ,q r2 ;X,Y n ) - X(q r2 ; Y n \X) 
-X(q r2 ;X) 

— T{q.ri , <2V 2 ! X, Y n ) — X(qv 2 ; X) 

= X(q ri ,q r2 ;X)+X(q ri ,q r2 ;Y n \X) 
-X{q r2 ;X) 

= X(q ri ,q r2 ;X)+X(q ri ;Y ri \X) 
+ l(qr 2 ; K n |g ri , X) - X(q r2 ; X) 
= Z(Qri , Qr 2 ;X) +X(g ri ; K n |X) -X(q r2 ;X) 

X)+^( ?n |X) 
-^( ?ri |F ri ,X)--H(g r2 )+^(g r2 |X) 
= X(g ri , g r2 ; X) +U (q ri |X) 



-H(g ri |r n )-H( ?r2 )+H(g r2 |X) 

D 



ln(l + a CF P s )+ln ( 1 + -i 



1 9 2 + Dj 



= In (1 + a CF P s ) 

/ (g 1 + P)(g 2 + P) \ 
+ 02(1 + 0^,)),/ 



= ln(l + a OT P s) + ln^^ rT -^ r j . (79) 

(a) and (d) follow from the fact that X i — > Y ri i — >■ g n is a Markov chain, and hence X(q ri ; X\Y ri , q r2 ) = 
and / H(q ri \Y n ,X) = H(q n \Y ri ). (b) and (c) follow from X(q r2 ; Y ri \X) = and X(q r2 ; Y ri \q ri , X) = 0, 
respectively, with respect to the Markov chain g r2 i — > X i — >■ y ri . ■ 
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With respect to the network symmetry and based on Lemma 111 one can express 

X(g r2 -Y r M=\n(l + a CF P s ) + \n ' ^ + D)( ^+ D > 



(80) 



D(l + a ri P s 

In order to have a successful transmission, the destination must first decode the relays signals and then 
X. From Eqs. (1671) to (17 lh . (1781) and (|80l) . to decode the relays signals at the detination, the following 
inequalities must be satisfied. 

h+D) (6 2 + D) 



ln(l + a CF P s ) + ln 
ln(l + a CF P s )+ln 



D(l + a r2 P s ) 
h + D) (6 2 + D) 
D{l+a ri P s ) 



<i? r <ln(l + aiP r ) , 
<R r <\n(l + a 2 P r ) , 



In 



[a ri +a r2 )P s + l 
D 2 



<2R r <ln((ai + a 2 )P r + l). 



(81) 



Therefore, the probability of decoding the relays signals at the destination is expressed as follows, 



Vn = Pr < max < In 



In (1 + a C jrP s )+ln 



\f{a ri + «r 2 ) P s + 1 
D 

'fix + D) (gg + D) 
D(l + a r _P s ) 



< Rr < 



min | In (^J {a[ + a 2 ) P r + lj , In (a min P r + 1)| L 



where a rmin = min {a n , a r2 } and a min = min {a x , a 2 }. 

After decoding the relays signals at the destination, the source signal is decoded subject to 

R < l(q ri ,q r2 ; X) = In (1 + a CF P s ) , 

where R — In (1 + P s s) is the source transmission rate. 
To summarize, we have shown the following. 

Theorem 6 The maximum throughput in the proposed CF scheme is expressed by 

^c,s= maxP c F acF ( S )ln(l + P sS ), 

S,D,Rr 



where acF and Vc are given by Eqs. (17 3l) and (|82l) . respectively. 



Analogously, Eq. (1851) yields the maximum continuous-layer expected-rate in this scheme. 



TQ C = maxVc 



D,R r 



- + 



J a 



ds. 



The integration limits are the solutions to F acF (so) 
respectively. 



(82) 



(83) 



(84) 



(85) 



s (1 + P s s ) / acF (so) and F acF (si) = s 1 / aoj? (si), 
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It turns out from the numerical results that the proposed CF scheme outperforms DAF and consequently, 
DF and AF, when the relay power to the source power ratio is higher than a threshold. 



Remark 3 If P s — > oo, Eq. (|73h is simplified to clcf ~ °' "nm+i) • ¥ Pr — > oo, then Vc — > 1 and 
acF ~ Q ri + a ro . In this high SNR asymptote at the relays, Eq. (85) meets the cutset-bound of Proposition^ 



in Section 



VII-A 



, and is optimum. 

VII. Upper-Bounds 

A. Cutset Bound 

The network cutset bound is the minimum of the maximum throughput and maximum expected-rate 
of the first-hop and the second-hop which lend itself to a closed form expression. The first-hop cutset is 
equivalent to a point-to-point single-input multiple-output (SIMO) channel with two receive antennas. The 
second-hop is equivalent to a multiple-input single-output (MISO) channel with two transmit antennas. 
The throughput cutset bound is the minimum of the maximum throughput in these two cutsets, that is 

^ Ss = maxe~ s (l + s)ln(l + Ps), (86) 

where P — min {P Sl P r }. 

Similarly, the maximum expected-rate of the diamond channel is upper-bounded by the minimum of 
the maximum expected-rates of those two cutsets, which is summarized below. 

Proposition 3 In the diamond channel, the cutset bound on the maximum expected-rate is specified by 

^ SiC =3E 1 (s ) -3E!( Sl ) - (s - l)e- S0 + (s, - l) e - Sl , (87) 
where s 1 = ^fi, and s = \J ^ A 2 - B 3 Ta + J - ^ with P = mm{P s ,P r }, A = 

1 1 1 _„ J D _1 I 1 

2P 6P 2 27P 3 ' U lU 3P ^ 9P 2- 

Proof- 
As showed in lfT5l . the maximum continuous-layer expected-rate is given by 

f°°— -sl'(s) 
11™ = max / Fjs)- ^-ds. (88) 

Noting Fg,(s) — (1 + s) e~ s based on EBl . we have 

7( S ) J 1 + Sl(s) 
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i . 



X Yr X r Y 



Fig. 3. Upper-bound model. 



The optimization solution to Eq. (189l) with respect to I(s) under the total power constraint P = min {P s , P r } 
is found using variation methods ll27l . By solving the corresponding Euler equation [|27l . we come up 
with the final solution as follows, 

7C = I 1 e~ s (1 +*)(-- 1 J ds, (90) 



where boundaries s and Si are the solutions to Psg + Sq — s — 1 = and s\ — s x — 1 = 0, respectively. 



Therefore, s x = ^fi, and s = s/vM 2 - £ 3 + A + - ^ with A = ^ - ^ - ^ , and 

5 = + . The indefinite integral (antiderivative) of Eq. (I90|) is 

e- s (1 + a) f - - 1 ) ds = (s - l)e~ s - 3^i(s). (91) 



s 

Applying the integration limits completes the proof. 



B. Relay-Cooperation (RC) Bound 

Here, a tighter upper-bound based on a full-cooperation between the relays is proposed. Let us define 
an upper-bound model by considering a full cooperation and power cooperation between the relays in the 
problem of interest (see Fig. |3|). The upper-bound model is equivalent to a dual-hop single-relay channel 
with two antennas at the relay. The following presents the throughput of this upper-bound model. 

Proposition 4 the maximum throughput in the above upper-bound model is given by 

P 
'P 



Ps 



n% c =m&x{l + s)[l + s-^ )e"H 1+ ^Jln(l + P s s) . (92) 



Proof: The optimum relaying strategy for dual-hop single-relay channels is DF. In the same general 
outline to the proof of Theorem [l[ let F (s) denote 



where a = -p-hQhJ and h 



F(s) = Pr{a r > s}Pr ja > psj , (93) 
hi hi ■ The maximum throughput in general can be expressed as 

n RC>s = maxF(s) ln(l + P s s). (94) 

s 

Analogously to the proof of Theorem [l[ we can restrict our attention to p = or p = 1, where p is the 
correlation coefficient between the signals transmitted from two relay antennas. To prove by contradiction, 
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first we assume that p° — 1; next we shall show that F p=0 (s°) > F p=1 (s°), which implies a contradiction 
and concludes p° = 0. Defining 

A 



g(s) = \n(l + P s s) — 
, , a F P =i{s) 



r s 



Eqs. (UJ) and (UJ) hold. 



f P =i(s) 



Noting 



we have 



F p=1 (s) = (l + s) e 



-si- 



rs 



1 + 5 



It can be shown that 



where 



( 1 + s )( 1 + &)- 1 



r(s) < s, Vs > s t , 



a VP, 2 + 4P s P r + 20P r 2 - P s + 2P r 



2P S + 4P r . 

Hence, Eqs. (Il9l) and (|2Q|) still hold by redefining r(s) and s t as above. 

As TZ(s) is a continuous function, one can conclude that < s° < s t . Noting 



s t < s c 



2W_i 



-1 



and 



P 



P p=0 (s) > F p= i(s), Vs < s c 



+ 1 1 § w 2.5129^ 



P s 



(95) 
(96) 



(97) 



(98) 



(99) 



(100) 



(101) 



(102) 



yields F p=0 (s°)> P p=1 (s°) and thereby, p° = and a = a± + a 2 . Substituting the channel gain CDFs in 
Eqs. (|93l) and (|94l) . the maximum throughput of the DF diamond channel is given by Eq. (1921) . 



The highest expected-rate of dual-hop single-relay channels has been studied in [17]. Here, only the 
final solution is mentioned as 



TZ^ = max / / f ar (t)F a (s) 

Is(a r ) J J 
l r (a\a r )0 

The power constraints at the transmitter and the relay are 



— sl' r (s\a r = 


t)ds 


1 + sl r (s\a r 


= t) 


= t) = P r . 





dt. 



(103) 



(104) 
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As the maximum transmission rate of the relay can not exceed its successfully decoded rate, the constraint 
on rate is 

sl' r {s\a r = t)ds f* sl' s (s)ds 
" /o 1 + sl s (s) ' 







$ If ( $ Qii" 



t) 



Vt. 



(105) 



The optimization problem of Eq. (11031) can be solved numerically using the algorithm proposed in [1301 . 

Following a similar outline in the proof of Theorem 2] and Proposition 4, one can show that the optimum 
transmission strategy at the relay is to transmit uncorrected equal power signal s from both of the relay 
antennas at each layer. Thus, F a _ r (s) = F a (s) = (1 + s)e~ s . Substituting in Eq. (11031) . we come up with 
the upper-bound as follows, which does not lend itself to a closed form formulation. 

Proposition 5 In parallel relay networks, the maximum expected-rate at the destination is bounded by 

-dsdt, (106) 



1Z r c c = max Ite 

I 3 (a r ) 
I r (a\a r )0 



00 00 

t , s(s + l)e~ s I' r (s\a r = t) 



1 + sl r (s\a r = t) 



subject to the power and rate constraints Eqs. (|104h and (|105r) . respectively. 



C. DF -Upper-Bounds 

As pointed out earlier, the continuous-layer coding for DF relaying can not be directly solved by 
variations methods. Here, two upper-bounds for the maximum continuous-layer expected-rate in DF 
scheme are obtained. Let us define a DF -upper-bound model as a diamond channel with uninformed 
transmitters, wherein the channel gains of the transmitter-relay links are both max{a n ,a r2 }, and that 
of the relay-destination links are and a 2 . This channel can be modeled by a dual-hop single-relay 
channel with the channel gains a r = max{a ri , a r2 ] and a for the transmitter-relay link and the relay- 
destination link, respectively. Clearly, the maximum expected-rate of this model yields an upper-bound 
on the maximum expected-rate of DF relaying. 

T he op timum relaying strateg y in th e DF -upper-bound model is decode-forward, and is given by 



Eq. (11031) . Analogous to Section 



VII-B 



, it can be shown that the optimum transmission strategy at the 



relay is to transmit uncorrected equal power signals from the relays at each layer. Hence, substituting 
F ar (s) = e _s (2 — e~ s ) and F a (s) = (1 + s)e~ s in Eq. (|103r) . we come up with the upper-bound as 
follows, which does not lend itself to a closed form formulation. 

Proposition 6 In the DF diamond channel, the maximum expected-rate at the destination is bounded by 



K 3 =,nax2^'(e- - l] fet 1)e "' J f ( ' k < 7 
is(a r ) J v 'J 1 + sl r (s\a r = t) 

Ma|a r )0 ^-0- 



dsdt, 



(107) 



subject to the power and rate constraints Eqs. (11041) and (11051) . respectively. 
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The cutset bound of the DF -upper-bound model results in a closed form expression. The results are 
summarized below. 

Proposition 7 The cutset bound of the DF -upper-bound model is specified by min {72.i,7^2}» where 

fti=4Ei (si) - 2Ei (2si) + e" 2si - 3e~ Sl - In (l- e~ Sl ) -0.1157, 

^ 2 =3E!(s 2 ) - (s 2 - l)e~ S2 - 0.1296. (108) 
si is the solution to — 2 ,~ e - s = 1 + P s Si and s 2 = V "\A4 2 — -B 3 + ^4 + ; ; = = — ts - w*^* 

4 _ 1 1 1 /7f7/7 R 1 I l 

~~ 2P r 6P r 2 27P3> " " 3P r 9P r 2- 

Proof: The bound on the second hop, i.e., 7?.2, is a direct result of Proposition ^ 

Noting F ar (s) = e~ s (2 — e~ s ) in the first hop, analogous to the proof of Proposition we have 

f°° -e- s (2 - e- s ) F(s) 

H x = max / ^ -. ( 1 } ds. (109) 

Ks) Jo 1 + sl(s) 



The optimization solution to Eq. (1109b with respect to J(s) under the total power constraint P s is found 



by solving the associated Eiiler equation 11271 . which leads to 



f Sl /2 1e~ s — 1 \ 

K l = / e-*(2-e-°)(- + T — zr )ds, (110) 



where boundaries s and si are the solutions to 



so 

2-e- s _ -. , j-) „„ , 2-e" s l 



1 + P s s and - , e _ s, = 1, respectively. 



o l-e" s a u 25i l-e" 3 l 



The indefinite integral (antiderivative) of Eq. (|110l) is 

- 4Ei(si) + 2Ei(2si) -e- 2si + 3e~ Sl +ln (l - e" Sl ) . (Ill) 
Applying the integration limits completes the proof. 



VIII. Numerical Results 



The achievable throughput, two-layer expected-rate, and continuous-layer expected-rate in the proposed 
multi-layer relaying schemes and their upper-bounds are shown respectively in Figs. I4] to H for P s = 
dB and < P r < 60 dB. Note that the rates are expressed in nats. When |f , namely powers ratio, is less 
than 25 dB, DAF is the best scheme. In higher values of the powers ratio, CF is the superior. AF has the 
worst performance for |f > 10 dB. When the powers ratio goes to infinity, CF meets the upper-bounds. 




Fig. 5. 



Two-layer expected-rate in the diamond channel. 
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P s = dB 

0.7 | 1 1 1 



0.6 - 




P T (dB) 

Fig. 6. Continuous-layer expected-rate in the diamond channel. 

IX. Conclusion 

The main goal of the paper is to propose simple, efficient, and practical relaying schemes to increase 
the average achievable rate at the destination in dual-hop parallel relay networks with Rayleigh block 
fading links and uninformed transmitters. To this end, different relaying schemes, in conjunction with 
the broadcast approach, were proposed. The performance of the proposed schemes were derived and 
numerically compared with two obtained upper-bounds. 
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